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Abstract 

We consider for the full time-dependent Maxwell's equations the inverse problem of 
identifying locations and certain properties of small electromagnetic inhomogeneities 
in a homogeneous background medium from dynamic boundary measurements on the 
boundary for a finite time interval. 
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1 Introduction 

I The ultimate objective of the work described in this paper is to determine locations 

and certain properties of the shapes of small electromagnetic inhomogeneities in a ho- 
mogeneous background medium from dynamic boundary measurements on part of the 
boundary and for finite interval in time. Using as weights particular background solu- 
' tions constructed by a geometrical control method we develop an asymptotic method 

based on appropriate averaging of the partial dynamic boundary measurements. 
For stationary Maxwell's equations it has been known that the Dirichlet to Neumann 
map uniquely determines (smooth) isotropic electromagnetic parameters, see [T6], jl8j. 
|20j . We will provide in this paper a rigorous derivation of the inverse Fourier transform 
of a linear combination of derivatives of point masses, located at the positions Zj of the 
inhomogeneities, as the leading order term of an appropriate averaging of (partial) dy- 
namic boundary measurements of the tangential components of electric fields on part 
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of the boundary. Our formulas may be used to determine properties (location, relative 
size ) of the small inhomogeneities in case a single, or a few (tangential) boundary 
electric fields are known. Our approach differs from [1], [2], [3], [1], [22] and is expected 
to lead to very effective computational identification algorithms. 
Our main result is given by: 

Theorem 4.1 Letrj e R*^. Let Ea be the unique solution in C^{0,T; X{Q))nC^ {0,T; L'^{n)) 
to the Maxwell's equations ^ with ip{x) = -q^e^^'^, ip^x) = —i^JjL^\ri\'q-^e'^^'^, and 
f{x,t) = j^-'-e«'?-2^-«Vw)l'?l*. Suppose that T and T geometrically control Q, then we have 

J J \_^^ ' ^ Ea X n — curl E x n) + dtO^^ ■ dt{ curl Ea x n — curl £' x n) da{x)dt ■ 

m 

a" ^(^0 - /ii)e2*''-^^M,(r?) • + ©(a^), 

where Qn is the unique solution to the Volterra equation /i20\) with defined as the 
boundary control in U8\) and Mj is the polarization tensor of Bj, defined by 

{Mj)k,i = efc • ( / {vj + (^ - l)^\ + {y))y ■ ei ds,{y)). 

Here (61,62) is an orthonormal basis of . The term 0{a^) is independent of the 
points {zj, j = 1, • • • , m}. 

For discussions on closely related (stationary) identification problems we refer the 
reader to [IS], [21], [6], and [TO]. 



2 Problem formulation 

Let Q be a bounded C^-domain in W^, d = 2,3. Assume that Q contains a finite 
number of inhomogeneities, each of the form Zj + aBj, where Bj C M"^ is a bounded, 
smooth domain containing the origin. The total collection of inhomogeneities is Ba = 
L)'jLi{zj + aBj). The points zj £ il., j = 1, . . . , m, which determine the location of the 
inhomogeneities, are assumed to satisfy the following inequalities: 

\zj — zi\ > Co > 0, V j / / and dist{zj , d^l) > cq > 0, V j. (1) 

Assume that q > 0, the common order of magnitude of the diameters of the inhomo- 
geneities, is sufficiently small, that these inhomogeneities are disjoint, and that their 
distance to M*^ \ is larger than co/2. Let and Eq denote the permeability and 
the permittivity of the background medium, and assume that hq > and Eq > are 
positive constants. Let > and ej > denote the permeability and the permittivity 
of the j-th inhomogeneity, Zj + aBj, these are also assumed to be positive constants. 
Introduce the piecewise-constant magnetic permeability 

[ fj-j, X £ Zj + aBj, J = 1 . . . m. 

If we allow the degenerate case a = 0, then the function no{x) equals the constant //q- 
The electric permittivity is defined by ea{x) = Eq, for all x £ Let n = n(x) denote 
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the outward unit normal vector to at a point on d^l, dtu = — and A means the 

ot 



Laplace operator defined by Au = 



dx 
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In this paper, we will denote by bold letters the functional spaces for the vector 
fields. Thus H^{fl) denotes the usual Sobolev space on Q and H'^(fi) denotes {H^{Q))'^ 
and L^(r2) denotes {L'^{Q))'^. As usual for Maxwell equations, we need spaces of fields 
with square integrable curls: 

H( curl ■,Q)={ue L2(S1), curl u E L^{n)}, 

and with square integrable divergences 

H(div ■,n) = {uG L'^in), div u E L'^{n)}. 

We will also need the following functional spaces: 

Y{n) = {u £ L'^{n),dw u = inn}, x{n) = ii\n) r\Y{n), 

and TL'^(dn) the space of vector fields on that lie in L^(9ri). Finally, the "minimal" 
choice for the electric variational space would be 

XAr(O) = E H( curl nH(div t; x n = ondn}. 

Now, we introduce the following time-dependent Maxwell equations (associated to 
the electric field) 



{Sad^ + curl — curl )Ea 



(3) 



where Ea E R'^ 
the initial data 



in n X (0,r), 

diY{eaEc) = mnx{0,T), 
Ea\t=o = dtEa\t=o = tp in ft, 

is the electric field, / the boundary condition for E^ x n, and ip and ■0 



Let E be the solution of the Maxwell's equations in the homogeneous domain: 

{eod'^ + curl — curl )E = in 17 x (0, T), 

div (eoE) = in 17 x (0,r), 
E\t=o = 93, dtE\t=o = ip in 



(4) 



£^ X n| 



ianx(o,T) 

Here T > is a final observation time and <f,ipe C°°(Q) and / E C'^{Q,T]C°°{dVt)) 
are subject to the compatibility conditions 

dff\t=o = (AV) X n\9n and +V|t=o = (A^) x n|an, ^ = 1, 2, . . . 

it follows that ^ has a unique solution E E C°°([0,T] x 17). It is also known (see for 
example jl7j ) that since 17 is smooth {C^— regularity would be sufficient) the non ho- 
mogeneous Maxwell's equations ([3]) have a unique weak solution Ea E C°(0, T; X(17)) n 
Ci(0,r;L2(17)). Indeed, curl ^« belongs to C°(0, T; X(17)) n Ci(0, T; L2(17)). 
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3 Asymptotic formula 

We start the derivation of the asymptotic formula for curl Ea x n with the following 
estimate. 

Lemma 3.1 The following estimate as a ^ holds: 

\\dt{Ea - -E')||ioo(o,T;L2(n)) + \\Ea- ^1 1 L°° (o,r;X]v(n)) < Co, (5) 

where the constant C is independent of a and the set of points {zj}^^ provided that 
assumption ([ip holds. 

Proof. From ([3])-([l]), it is obvious that Ea — EG X]y{i^), then due to the Green 
formula we have for any v G X]\f{^}): 

I eod^{Ea -E)-^ dx+ [ — curl {E^ - E) ■ curl v dx = (6) 

m ^ 1 /■ 

> ( ) / curl E ■ curl v dx. 



Let Vq be defined by 



v« G Xn{^), 

curl curl Vq, = dt{Ea — E) in Vl. 



Then, 

/ — curl {Ea - E) ■ curl ^ra dx = - [ dt{Ea - E) ■ (S„ - E) dx 
Jn fJ-a Jn 



-Idt I \Ea-E\^ dx 



n 



2 

and by Green formula, relation ([7|) gives: 



/ df{Ea — E) ■ dx = [ curl — curl dt^a ■ dx 
Jn Jn l^a 



— curl ^t^ra ■ curl dx 
n /^a 

1 f 1 

-dt I — I curlv^p dx. 

2 Jn fJ-a 



Thus, it follows from ([6]) that 

1 



^odt I — I curl dx + dt / \Ea — E]"^ dx = 
Jn f^a Jn 

™ 1 If 
2N ( ) / curl ii^ • curl Va dx. 

7r[ /^O fJ-j Jzj+aBi 



(7) 
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Next, 



m ^ 1 /■ 

I V'( ) / curlS- curl Vol < C|| curl £'||L2(f5^)|| curl Vo||l2(q). 

Since E G C°°([0,r] x H) we have 

m 

II curlE||L2(g^) < II curlE||ioo(g^)a(^|Bj|)2 < Ca, 



which gives 



" 1 1 /■ 

> ( ) / curl • curl Vfjcij; I < Ca 1 1 curl Va 1 1 l2 



and so, 



eo^t / —I curl v„p dx+(9t / \Ea-E\^ dx < Ca{ [ —\ curl v„|^ dx+ [ |Sa-£^P dx)^/^. 

(8) 

From the Gronwall Lemma it follows that 

{[ —\ curl V«|2 ^ ( /■ _ ^^)l/2 < (g^ 

Combining this last estimate ([9|) with the fact that 

\\dtiEa - -E)||^oo(o,T;H-i(f^)) - ^11 C^l^l '^alli°°(0,T;L2(C)) 

the following estimate holds 

ll-E^a - -E^o||l°°{0,T;L2{C)) + \ \dtiEa - -E^o) I |l°° (0,T;L2(n)) < Co. (10) 

Now, taking (formally) v = dt{Ea — E) in ^ we arrive at 

eodt I [\dt{E^ - E)\^ + —I curl {E^ - E)\^] dx = 
Jn /^a 

m 

2y ( ) / curl E ■ curl dt{Ea, - E) dx. 

By using the regularity of ii^ in 17 and estimate (jlOp given above, we see that 



y.{—-—)! cmlE- cnildtiEo^-E) dx\ < C| | curl E| |h2(b.)| |5f (E„ - E)| |h 



where C is independent of t and q, and so, we obtain 



dt I [\dt{Ec, - ^)P + —I curl {E^ - E)\'^] dx < Ca^ 
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□ 



which yields the fohowing estimate 

\\dt{Ea - -E')||ioo(o,T;L2(f7)) + " {0,T;Xn{^)) - ^'o, 

where C is independent of a and the points {zj}JLi. 

Now, we can estimate curl — curl Eq as follows. 

Proposition 3.1 Let E^ and E he solutions to the problems and ^ respectively. 
There exist constants < ao, C such that for < a < i^e following estimate holds: 

II curl {Ea - £^o)||Lo°{0,T;L2(n)) < Ca, (11) 
Proof. To prove estimate (|lip it is useful to introduce the following function 

j-T 

v{x)= v{x,t)z{t) dt e L'^in), (12) 
Jo 

where v € L^{0,T;L'^{n)) and z{t) is a given function in Cq°(]0,T[). 
Then, 

E{x) = [ E{x,t)z{t) dt and Ea{x) = [ Eaix,t)z{t) dt £ X{n), 
Jo Jo 

which by relation ([5]) give 

f {Ea-E)£ilHn), 

curl curl (E^ — E) = 0(a) in 
div {Ea - E) = inn, 
, {Ea - E) X nlen = 0, 
and so, 

1 1 curl (£;„-£;) 1 1 L2 (f,) =0 (a). (13) 

The fact that curl {E^ - E) belongs to L°°(0, T; V'{n)) and by using estimate ^ we 
deduce that 

/ \ cnil Ea{x,t) - cnT\E{x, t)\'^ dx = 0{a^) a.e. in t G (0, T), 
Jo. 

which means that 

II curl [Ea — -E')||L2(r2) = 0{a) a.e. in t G (0,T). 

Thus, estimate (llip follows immediately if we take the sup on t G (0, T) in the last 
relation. □ 
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Before formulating our main result in this section, let us denote 
the unique vector-valued solution to 



1, . . . ,m 



in Bj , and 



<I>o is continuous across dBj, 



(14) 



lim \^j[y)\ 

( — >+oo 



0, 



where Vj denotes the outward unit normal to dBj, and superscripts — and + indicate 
the limiting values as the point approaches dBj from outside Bj, and from inside Bj, 
respectively. The existence and uniqueness of this <I>j can be established using single 
layer potentials with suitably chosen densities, see |6J for the case of conductivity 
problem. For each inhomogeneity Zj + aBj we introduce the polarizability tensor 
Mj which is a d X d, symmetric, positive definite matrix associated with the j-th 
inhomogeneity, given by 



{Mj)k,i = efc • ( / {vj + { 



dB-i 



>0 



iy))y • ei d(Tj{y)). 



(15) 



Here (ei, . . . , e^) is an orthonormal basis of W^. In terms of this function we are able 
to prove the following result about the asymptotic behavior of curl Ea ■ i^j\d{zj+aBj)+- 



Theorem 3.1 Suppose that ([I]) is satisfied and let ^j,j = l,...,mbe given as in ((j^[ ). 

Then, for the solutions E^, E of problems and respectively, and for y € dBj 
we have 

( curl Ec,{zj + ay) ■ yj)\d{z,+aB,)+ = curl E{zj,t) ■ Vj (16) 
+(l-^)^| + (y)- curl i?(z„t) + 0(1). 

/iO OVj 

The term o(l) uniform in y £ dBj and t G (0, T) and depends on the shape of {Bj}JLi 
and Q, the constants cq, T, /xq, {fJ-jYjLi, the data ip,ip, and f, but is otherwise inde- 
pendent of the points {zj}'jLi. 

Proof. 

Let £a = curl Ea{x,t) and Sq = curl E{x,t). Then, according to ©-(jl]) we have 

eod^Ea — curl —£a = and curl £a = 0, for x e 0. (17) 

We restrict, for simplicity, our attention to the case of a single inhomogeneity, i.e., 
the case m = 1. The proof for any fixed number m of well separated inhomogeneities 
follows by iteration of the argument that we will present for the case m = 1. In order 
to further simplify notation, we assume that the single inhomogeneity has the form 
aB, that is, we assume it is centered at the origin. We denote the electromagnetic 
permeability inside aB by and define the same as ^j, defined in (I14j) . but with 
Bj and fij replaced by B and fi^, respectively. Define v to be the outward unit normal 
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to dB. Now, following a common practice in multiscale expansions we introduce the 

X ~ Q 

local variable y = then the domain Q = (— ) is well defined. 

a a 
Next, let w be given in C^{]0,T[). For any function v e L1(0, T; L2(0)), we define 



v[x) 



T 



v{x,t)zu{t) dt e L^(0). 



We remark that dtv{x) 
satisfies 



v{x,t)w' (t) dt. So that we deduce from (fTTl) that £a 



1 



curl — £a 



Ea vj"{t) dt in Q, 
curl £n = in 17. 



Analogously, f satisfies 



curl<S = / Ew"{t) dt in 



1 

/^O JO 

curl £ = in 0. 



Indeed, we have £o,xn = £xn = curl qqJ x n on the boundary 50, where curl qq 
is the tangential curl. Following [4j and [Ij, we introduce (7* as the unique solution to 
the following problem 

f - Q _ 

Ag*=0 inJ7 = (-)\Sandin5, 
a 



is continuous across 

= -{^iQ- H^)£{ay) ■ V ondB, 



dql dql 



U Of 

q*^ = on 017. 



The jump condition 



Mo- 



du 



du 



-{pLQ — ijL^)£{ay) ■ u on dB 



guarantees that £a{x) — £{x) — grad ^^^(f ) belongs to the functional space Xi\]{Q), 
where grad is the tangential gradient. Since 



curl — {£a-£- grad g*(-)) 

IJLa "a 



Ea - xi^ \ aB)E + —x{aB)E] zu"{t) dt in O, 



curl {£a — £ — grad yQai—)) = in 17, 



a 



^ {£a-£- gradyg*(f)) X n = 



on 



on. 
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where xi'-^) is the characteristic function of the domain uj, we arrive, as a consequence 
of the energy estimate given by Lemma l3.lt at the following 

(4-^- grad^<?*(f))EX^(l^), 

1 - - X 
curl — (£a — £ — curl vq*( — )) = 0(a) in 17, 
fia a 

curl {£a- £ - grad yqai — )) =0 in 
{£a-£ - grad ^g* (f )) x n = on dQ. 
From ^ we know that this yields the following estimate 

II curl — {£a -£ - grad yq*a{-))\\L^(n) + ll-^a - - grad yq*ai-)\\L^n) < Ca, 



and so. 



{£a-£ - grad „g* (-)) • = 0(a) on d{aB). 



a 



Now, we denote by q^, be the unique (scalar) solution to 

'Ag* = inM'^\BandinS, 
g=K is continuous across dB^ 



Mo 



dq* , 
5z/ 



9zv 



-{fiQ — fi*)£{0) ■ u ondB, 



lim = 0. 

. |yH+oo 



In the spirit of Theorem 1 in [6j it follows that 



grad q^ - grad „g* )(-)| Ilz^) < Ca^/^ 



which yields 



{£a — £ — grad yq*{—)) ■ v = o(l) on d{aB). 



a 



Writing q^, in terms of <I>* gives 
T 

L 



curlS„(ay)-i/)|a(„B)+-z.- curl ^(0, t)-(^-l)^| + (y)- curl ^(0, t)] ro(t) = o(l). 



for any G C[^(]0,T[), and so, by iterating the same argument for the case of m 
(well separated) inhomogeneities Zj + aBj,j = l,...,m, we arrive at the promised 
asymptotic formula ([16]) . 



□ 



4 The identification procedure 

Before describing our identification procedure, let us introduce the following cutoff 
function P{x) S C^{i}) such that /? = 1 in a subdomain of that contains the inho- 
mogeneities Ba and let rj G M.^. We will take in what follows E{x,t) = 7y-Lg«f?-a:-jVw)l'?l* 
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where rj-^ is a unit vector that is orthogonal to rj which corresponds to taking (p{x) 



-i^/JIo\ri\r] e^"^'^ , and f{x,t) = r] x 



ngW^ ^VJ^lvlt and assume that we 



r] e 

are in possession of the measurements of: 

curl EaXn on T x (0,r), 

where F is an open part of dQ. Suppose now that T and the part T of the boundary 
are such that they geometrically control which roughly means that every geometrical 
optic ray, starting at any point x E at time t = hits T before time T at a non 
diffractive point, see [5j. It follows from [17] (see also [l3], [H] and [l2]) that there exists 
(a unique) gr/ € -f^o (0, T; rL^(r)) (constructed by the Hilbert Uniqueness Method) such 
that the unique weak solution to 

{d^ + curl curl )wr, = in x (0, T), 
div'u;^ = inl7x(0,r), 

Wr,\t=o = Pix)v-^e'''-'',dtWr,\t=o = inn, (18) 

Wr, X n|rx(0,T) = 9r], 

satisfies Wr^{T) = dtWrjiT) = in 17. 

Let 6*^ G i7i(0,r;TL2(r)) denote the unique solution of the Volterra equation of 
second kind 



dter,{x,t) + 



-i\r,\{s-t) 



{Offix, s) — i\r]\dtOrj{x, s)) ds = g^ix, t) for x G F, t G (0, T), 



9^(x,0)=0 forxGF. 

(19) 

The existence and uniqueness of this 9^^ in if^(0, T; rL^(r)) for any rj G M'^ can be 
established using the resolvent kernel. However, observing from differentiation of (jl9p 
with respect to t that O-q is the unique solution of the ODE: 



= S'^\*^t{e-'\'^\^gr^) for x G F, t G (0, T), 



6l^(x,0) = 0,9i6l^(x,r) = forxGF, 



(20) 



the function 6^ may be find (in practice) explicitly with variation of parameters and it 
also immediately follows from this observation that On belongs to H'^{0,T;TL^(T)). 
We introduce as the unique weak solution (obtained by transposition as done in 
|15j and in [II] [Theorem 4.2, page 46] for the scalar function) in C^{0,T; X{Q,)) n 
C^{0,T;L'^{n)) to the following problem 

[d^ + curl curl )u^ = in 17 x (0, T), 



div Vrf 
Vrj\t=0 



infix (0,r), 
in n, 

m 



.^B.) G Y(Sl) in 0, 



V-q X n|gQx{0,T) = 0. 
Then, the following holds. 
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Proposition 4.1 Suppose that T and T geometrically control Q. For any t] £ M we 

have 

r [ Qrj-i curl Vr, X n) da{x)dt = a2V/io(l - ^)e^''T'^r] ■ [ {vj (21) 
h Jr A^o JdBj 

+(^ - l)^\ + {y))v ■ y ds,{y) + o(a2). 

Proof. Multiply the equation {d^ + curl curl )vrf = hy Wr^ and integrating by parts 
over (0, T) X Q, for any r/ E M"^ we have 



m „ 
~^ ^0 JdBj /^O 



^ / 9ri ' { curl X n) da{x)dt. 
Jo Jr 



Now, we take the Taylor expansion of ae*"^'^' in the left side of the last equation, we 
obtain the convenient asymptotic formula (j2ip . 

□ 

To identify the locations and certain properties of the small inhomogeneities Ba let 
us view the averaging of the boundary measurements 

curl Ea X n|rx(o,T), 

using the solution 0^ to the Volterra equation (fT9|) or equivalently the ODE ([20]) , as a 
function of rj. The following holds. 



Theorem 4.1 Letr]£R'^. Let Ea be the unique solution in C^{0,T; X{Q))nC^ {0,T; L'^{Q)) 
to the Maxwell's equations ^ with ip{x) = ry-'-e*'^'^, 'ip{x) = — iy^|ry|r/-'-e*'^'^, and 
f{x,t) = r/-'-e*''■^~*^/'^l''l^ Suppose that T and T geometrically control O, then we have 



J J \_^'n ' ( c^'^l Ea X n — curl E x n) + dtOr^ ■ dt{ curl Ea x n — curl £' x n) da{x)dt 

2 



(22) 

where On is the unique solution to the Volterra equation i20\) with g^ defined as the 
boundary control in U^) and Mj is the polarization tensor of Bj, defined by 

iM,)k,i = e,. • ( / + (^ - l)^| + (y))y . ei ds.iy)). (23) 

JdBj ^0 Ol/j 

Here (61,62) is an orthonormal basis ofM.'^. The term O(a^) is independent of the 
points {zj, J = 1, • • • , m}. 
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Proof. From dtO^iT) = and ( curl x n— curl E x n)|f=o = the term / / dtO^^ ■ 

Jo Jr 

dt{ curl Ea X n — curl E x da{x)dt has to be interpreted as follows 

ndtOr) ■ dt{ curl EaXn- curl E xn) = - / 8^9,^ ■ ( curl EaXn- curl E x n). 
Jo Jr 

(24) 

Next, introduce 



Ea^r,{x,t) = E{x,t) + I e-'^^'i^'vr,{x,t-s)ds,xen,te{0,T). (25) 



We have 



'0 Jr 

cT 




'0 Jr 

l-T 



■ ( curl EaX u.— curl £' x n) + dtOr, ■ dt{ curl E^ x n — curl £^ x n) 



( curl Ea X n — curl E^^r] x n) + dtOn • dt{ curl E^ x n — curl Ea,r) x n) 



+ 




Jr 



u^(x, t — s) X 



nds + dtOr, ■ dt [ e-*v^l^l"w^(x, t - s) x n ds 
Jo 



Since 9^^ satisfies the Volterra equation ([20]) and 



dt{ f e-'^\'^\'v^{x, t-s)xnds) = ^^(-e-^v^l^l* f e^v^l^l"t;^(x, s) x n ds) 
Jo Jo 

= i^|r?|e"^v^l^''l* / e*v^l''l"'t;^(x, s) x n ds + Vr^{x, t) x n, 
Jo 

we obtain by integrating by parts over (0, T) that 
rT r ^ rt 




'0 -^JT 
fT 







e -'^\'^\'v^{x,, 



t-s)xnds + dtOrj ■ dt / e^'^^'^^'vr^ix, t - s) x n ds 



Jr 



iv^ix,t) X n) • {dtO^ + / e^(s)e*v^l''l(*-") ds) 




-i^\r]\{e-'y'^"'^''^'dtO^{t)) ■ / e^^l''l^i;^(x,s) x n ds dt 

Jo 

rT 

Vr,ix, t)xn- (dtO^ + / iOr^is) - i^\v\dt9^is))e'^\''\^'''^ ds) dt 
lo Jr Jt 

— I I 9ri{x^t) ■ ( curl Vri{x,t) X n) dt 
Jo Jr 

and so, from Proposition 14.11 we obtain 

^6rj ■ ( curl EaXn— curl £^ x n) + dtOr^ ■ dt{ curl EaXn— curl E x n) 



Jr 



5^(1 - ^)e'^^-^r^ . I (., + (^ - l)^| + (y))^ • y dsM 
j=i ^0 -^es. 




Jr 



( curl Ea X n — curl E^^r] x n) + dt9n • dt{ curl Ea x n 



curl Ea^n x n) + o(a^) 
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In order to prove Theorem 14.11 it suffices then to show that 



rT 



10 JT 

Since 



9n-{ curl E^qXii— curl Ea^rjXn)+dtOrj-dt{ curl EaXn— curl i^Q^^xn) = o( 



a 



(26) 



(d? - curl — curl )( / e-'^^'^^'vJx, t - s) ds) 
/^o Jo 



^,(l_^)^x (i., + (^-l)^| + (y))e^^-^5e(,^+„5^)e-v^l^l* inOx (0,r), 



(/ e-*v^l^l%(x,t-s)(is)|t=o = 0,9t(/ e-*v^l^l%(x,t-s)ds)|t=o = in 17, 

JO JO 

( f e-'^^''^'vr,{x,t- s) ds) X n|9Ox(0,T) = 0, 



JO 

it follows from Theorem 13.11 that 

{d^ - curl — curl ){Ec, - S„ ^) = o(a^) in x (0, T), 

^0 

(E'q, — Ea,-q)\t=0 = 0,dt{Ea — Ea,n)\t=o = in 

(E'a - Ea,T]) X n|aQx(0,T) = 0- 

Following the proof of Proposition 13. H we immediately obtain 

ll^a - Ea,rj\\mn) = o(a^),te {0,T),x G Q, 

where o(a^) is independent of the points {zj}^-^^. To prove (j26p it suffices then from 
(j24p to show that the following estimate holds 



II curl £'„ X n - curl Ea,ri x n| |i2(o,r;TL2(r)) = o(a^). 
Let 6 be given in Cq° (]0,T[) and define 



E^,^ix)= / Ea,^{x,t)e{t) dt 



and 



Ea{x) = / Ea{x,t)e{t) dt. 

Jo 



From definition (1251) we can write 



{E^-E^)£-ll\Q), 
curl curl (E„ - = 0(q) G in 
div {E^ -Ea) = in 
{Ea - Ea) X n|an = 0. 



(27) 
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In the spirit of the standard elhptic regularity [9] we deduce for the boundary value 
problem ()27p that 

II curl {Ea - Ea) X n||L2(r) = O(a^), 
for all e e C^(]0,r[); whence 

1 1 curl {Ea — Ea) X n| |L2(r) = o(a^) a. e. in t E (0, T), 

and so, the desired estimate (|22p holds. The proof of Theorem 14.11 is then over. □ 

Our identification procedure is deeply based on Theorem 14.11 Let us neglect the 
asymptotically small remainder in the asymptotic formula ()22p . and define ^a{'n) by 

^a(^) = / / k • ( curl {E^ - X n) + 9*0^ • dt{ curl {E^ - E) x n) 

Jo JT L 

Recall that the function e^*'''^^' is exactly the Fourier Transform (up to a multiplicative 
constant) of the Dirac function 5-22^ (a point mass located at —2zj). From Theorem 
14.11 it follows that the function e^^"^'^^ is (approximately) the Fourier Transform of a 
linear combination of derivatives of point masses, or 

m 

i=i 

where Lj is a second order constant coefficient, differential operator whose coefficients 
depend on the polarization tensor Mj defined by (p3|) (see [6j for its properties) and 
^aiv) represents the inverse Fourier Transform of i^ai^)- The reader is referred to [6] 
for properties of the tensor polarization Mj. 

The method of reconstruction consists in sampling values of 'i^aiv) some discrete 
set of points and then calculating the corresponding discrete inverse Fourier Transform. 
After a rescaling the support of this discrete inverse Fourier Transform yields the loca- 
tion of the small inhomogeneities Ba ■ Once the locations are known we may calculate 
the polarization tensors (Mj)^-^ by solving an appropriate linear system arising from 
(I22|) . This procedure generalizes the approach developed in [3] for the two-dimensional 
(time-independent) inverse conductivity problem and generalize the results in [lj to 
the full time-dependent Maxwell's equations. 



5 Conclusion 

In this paper, we are convinced that the use of approximate formulae such as (122p 
represents a very promising approach to the dynamical identification of small inho- 
mogeneities that are embedded in a homogeneous medium. We also believe that our 
method yields a good approximation to small amplitude perturbations in the electro- 
magnetic parameters (for the example of electric permittivity ea(x) = eo + a£{x)) from 
the measurements: 

curli?„xn onrx(0,r). 

Our method may yield the Fourier transform of the amplitude perturbation e(x). This 
issue will be considered in a forthcoming work [7]. 
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